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Swift/XRT monitoring of the Supergiant Fast X-ray Transient 
IGR J18483-0311 for an entire orbital period 
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ABSTRACT 

IGR J18483— 0311 is an X-ray pulsar with transient X-ray activity, belonging to the new 
class of High Mass X-ray Binaries called Supergiant Fast X-ray Transients. This system is 
one of the two members of this class, together with IGR J 1 1215-5952, where both the orbital 
(18.52 d) and spin period (21 s) are known. We report on the first complete monitoring of 
the X-ray activity along an entire orbital period of a Supergiant Fast X-ray Transient. These 
Swift observations, lasting 28 days, cover more than one entire orbital phase consecutively. 
They are a unique data-set, which allows us to constrain the different mechanisms proposed 
to explain the nature of this new class of X-ray transients. We applied the new clumpy wind 
model for blue supergiants developed by Ducci et al. (2009), to the observed X-ray light 
curve. Assuming an eccentricity of e = 0.4, the X-ray emission from this source can be 
explained in terms of the accretion from a spherically symmetric clumpy wind, composed of 
clumps with different masses, ranging from 10 18 g to 5xl0 21 g. 

Key words: X-rays: binaries - X-rays: individual (IGR J18483— 031 1) 



1 INTRODUCTION 

The X-ray transient IGR J 18483— 031 1 was discovered during the 
observations of the Galactic plane with INTEGRAL in April 2003 
dChernvakova et al]l2003h . when it reached a flux of 10 mCrab in 
the 15^10 k eV energy ranee. Fiv e more hard X-ray outbursts were 
reported by iSguera et alj ( 12007b with INTEGRAL, three of which 
exceeded one day in duration. The strongest outburst (~1.8 days) 
occurred in April 2006 and reached a flux of 120 mCrab. The 
IGR 118483-03 11 broad band joint JEM-X/ISGRI INTEGRAL 
spectrum (3-50 keV) was fitted with an absorbed power law with 
a photon index, Y — 1.4 ± 0.3, a high absorption, Ah = 9^ X 
10 22 cm -2 (higher than the Galactic a t 1.4 x 10 22 cm -2 ), and a 
cutoff at - 22 keV dSguera et al]|2007h . 

A periodicity at 18. 55 ± 0.03 days was d iscovered in the 
ASM1RXTE light curves dlxvine & Corbetll2006t), and wa s inter - 
preted as the orbital period of a binary svstem. lSguera et all d2007h 
confirmed a similar period in the 20-40 keV IS GB1I INTEGRAL 
data (18.52 ± 0.01 days) and discovered pulsations at 21.052 6 ± 
0.0005 s with the X-ray monitor JEM-X. fslrinner et al.1 J2008h re- 



fined the orbital period value to 18.518 ± 0.005 d, by using the 
much denser sampling provided by the Swift/BKT data. 



From the observed values of th e orbital and pulse periods, 
the position in the Corbet diagram (Corbet 1986) suggested at 
first a Be/X-ray transient nature dSguera et al.ll2007h . but the op- 
tical and in frared observations of the X-ray error box estimated 
with Swift dSguera et al.1 120071) revealed that the donor star is a 
blue supergiant (B0.5 Ia), and not a Be star, located at a distance 
of 3-4 kpc (Rahoui & Chatv 2008). The X-ray position was later 
refined with Chandra dGiunta et alJ,l2009h . confirming this opti- 
cal/IR association. This implied the identification of this source as 
a new m ember of the class o f the Supergiant Fas t X-ray Transients 
(SFXTs; ISguera et al.1 [2001, ISguera et alj|2006l . iNegueruela etall 
l200d [Sidoli 20091) . although the dynamical range of its X-ray 
emission seems to be smaller than in other members of the same 
class. IGR J18483— 0311 and IGR 111215-5952 are the only two 
SFXTs where both orbital and pulse periods have been disc overed 
dSwank et alj200l[sidoli et alj2006llRomano et alj2009bl) . 
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Figure 1. Swift/XKT 0.2-10keV light curve of IGR J18483— 031 1 during our monitoring program, background-subtracted and corrected for pile-up, PSF 
losses, and vignetting. Downward-pointing arrows are 3<r upper limits. The upper limit with the wide symbol centered on MJD 55007.9 (~ 0.004 counts s —1 ) 
is obtained accumulating the 4 observations (individual upper limits) between MJD 55006.5 and 55010.2. Different colours mark different observations (see 
Table[T}, with a colour scheme that generally mimics the phase (top axis) with a P = 18.52 days iSeuera et alj2007l) . The inset zooms on observation 024. 



Table 1. Summary of the Swift/XKT observations. 



Seq. Start time (UT) End time (UT) Exp. Phase" 

(s) 



004 


2009 


06 


11 14:42:17 


2009 


06 


11 21:13:56 


1606 


0.07 


005 


2009 


06 


12 16:09:36 


2009 


06 


12 19:34:56 


1522 


0.12 


006 


2009 


06 


13 18:08:11 


2009 


06 


13 21:26:57 


1774 


0.18 


007 


2009 


06 


14 04:55:39 


2009 


06 


14 08:20:57 


1671 


0.21 


008 


2009 


06 


15 10:03:12 


2009 


06 


15 23:02:58 


1977 


0.28 


009 


2009 


06 


16 07:10:02 


2009 


06 


16 12:08:58 


1857 


0.32 


010 


2009 


06 


17 00:42:44 


2009 


06 


17 04:12:56 


2047 


0.36 


012 


2009 


06 


19 08:37:38 


2009 


06 


19 12:02:56 


2189 


0.49 


013 


2009 


06 


20 13:32:09 


2009 


06 


20 16:57:56 


2291 


0.55 


016 


2009 


06 


23 07:51:16 


2009 


06 


23 12:53:56 


2003 


0.70 


017 


2009 


06 


24 09:32:31 


2009 


06 


24 14:35:56 


2711 


0.76 


018 


2009 


06 


25 04:33:46 


2009 


06 


25 12:52:56 


1870 


0.81 


019 


2009 


06 


26 15:43:27 


2009 


06 


26 19:06:57 


1884 


0.88 


020 


2009 


06 


27 04:33:51 


2009 


06 


27 07:56:57 


1824 


0.91 


021 


2009 


06 


28 03:16:17 


2009 


06 


28 08:14:56 


1922 


0.96 


022 


2009 


06 


29 11:30:17 


2009 


06 


29 11:47:57 


1057 


0.03 


023 


2009 


06 


30 06:47:17 


2009 


06 


30 11:54:57 


2616 


0.08 


024 


2009 


07 


01 02:17:22 


2009 


07 


01 12:00:57 


2070 


0.13 


025 


2009 


07 


02 10:22:32 


2009 


07 


02 13:23:57 


2408 


0.19 


026 


2009 


07 


03 07:01:47 


2009 


07 


03 10:31:58 


2106 


0.24 


027 


2009 


07 


04 07:05:17 


2009 


07 


04 10:36:27 


2560 


0.29 


029 


2009 


07 


06 02:15:18 


2009 


07 


06 08:53:58 


1838 


0.39 


030 


2009 


07 


08 07:32:02 


2009 


07 


08 07:37:56 


340 


0.51 



a Calculated according to Sgue ra et alj j2007l) . 



2 OBSERVATIONS AND DATA REDUCTION 

The observations of IGR J18483— 031 1 were obtained as a Target 
of Opportunity (ToO) monitoring program with Swift. As shown in 



TableQ] the ToO observations started on 2009 June 1 1 with ~ 2 ks 
per day. The campaign lasted 28 days divided in 23 observations 
for a total on-source exposure of ~ 44 ks. 

The XRT data were processed with standard procedures 
(XRTPIPELINE vO. 12.1), filtering and screening criteria by using 
FTOOLS in the Heasoft package (v.6.6.1). Both WT and PC 
events were considered. The selection of event grades was 0-2 and 
0-12, for WT and PC data, respectively dBurrows et al]l2005h . We 
corrected for pile-up when required. The light curves were also cor- 
rected for PSF losses, vignetting and background-subtracted. An- 
cillary response files were generated with XRTMKARF, to account 
for different extraction regions, vignetting, and PSF corrections. 
We used the spectral redistribution matrices vOll in CALDB. 

The UVOT observed the target simultaneously with the XRT 
with the v filter (observations 004 through 006), and with the 
u filter (observations 007 though 010). For the remainder of the 
campaign IGR J 18483— 0311 was observed with the 'Filter of the 
Day', i.e. the filter chosen for all observations to be carried out 
during a specific da y in order to minimize the filter wheel usage 
dRoming et ai]|2005h . The data analysis was performed using the 
UVOTIMSUM and UVOTSOURCE tasks included in the FTOOLS. 
The latter task calculates the magnitude through aperture photom- 
etry within a circular region and applies specific corrections due 
to the detector characteristics. The reported magnitudes are o n 
the UVOT photometric system described in IPoole et all d2008l) . 
and are not corrected for Galactic extinction. At the position of 
IGR J18483— 031 1, no detection was achieved down to a limit of 
v > 21.07 mag and it > 21.19 mag. 

All quoted uncertainties are given at 90% confidence level 
for one interesting parameter unless otherwise stated. The spectral 
indices are parameterized as F v oc v~ a , where F v (erg cm~ 2 s^ 1 
Hz - ) is the flux density as a function of frequency v\ we adopt 
r = a + 1 as the photon index, N(E) oc E~ T (ph cm -2 s" 1 
keV" 1 ). 
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Table 2. XRT spectroscopy. 



Seq. Power-law Black body 

N H T Flux" X 2 Icd /dof N H kT R BB Flux" L b X 2 cd /dof 

(10 22 cm- 2 ) (2-10keV) C-stat (%) c (10 22 cm~ 2 ) (keV) (km) d (2-10keV) (2-10keV) C-stat (%) c 



005 




1 4+0-9 
1 -0.8 


1 90+ 037 

1,yu -0.94 


102.1(48.3) 


o c + 2.2 
d *°-1.8 


1 6+°' 6 


14+0- 03 
u,1 ^-0.01 


1 46+ 015 
1, ^ u -0.83 


0.2 


100.9(31 


90) 


006 


6.3±i;l 


i 7 +0.4 


15 46+0- 54 


1.1/34 


1+O.8 
°- x -0.7 


1 6+ ' 2 


43+0- 22 
U '^ J -0.14 


O+0.79 
°— 1 99 


1.2 


1.1/34 




007 


4 o+3.8 


1 n+0.9 

± -°-0.T 


9 7 9+0.39 
z - ' -0.92 


1.2/9 


9 c+2.0 
z.0_ 1 2 


1 6+°' 5 
i- D -o.3 


17+ 004 
u ' 1 '-0.02 


2 06+° :24 
z - UD -0.99 


0.2 


0.9/9 




010 


6 6+ 2 ' 7 

"■"-2.2 


n i +0.8 
—0.7 




119.0(27.2) 


3 5+ 1 ' 7 
■3-0-1.4 


1 q+0.3 
1,o -0.2 


n 10+0. 05 


1 oo +0 - 10 


0.1 


117.5(10 


58) 


012 


6 .4i» : ; 


1 s+o- 7 


2.78+8;" 


1.2/13 


3 2+ 17 


1 K+0.3 
i -°-0.2 


20+0' 05 
UiZU -0.02 


1 98+ ' 16 


0.2 


1.1/13 




022 


16.51^ 


2.3+1- 2 


5 79+0-33 
°-' a -4.33 


97.43(26.7) 


1 7+ 5 - 2 


1 c+0.6 
l-°-0.3 


r\ 97+O.28 
u,z '-0.08 


, 0+0.22 
o.<±o_2 24 


0.4 


98.02(12 


17) 


023 


6.6+1* 


I 4+0-3 
-0.3 




0.9/51 


7+O.8 
°-'-p.7 


1 7+0.2 
'-'-0.2 


u - oa -o.io 


13.26+ ,' 62 


1.4 


1.0/51 




024 


10.4+Jj 


i rr+0.5 
± -°-0.4 


19.97+^ 


0.7/25 


6 3+ 18 
— 1.4 


1 9+°- 3 
^-o.s 


07+0.19 


14.78+.°;?* 


1.6 


0.7/25 




025 


71 +i.i 
'•i-i.o 


1 6+°- 3 


09 OO+0.87 
oi ' OJ -4.92 


1.1/71 


o+0.7 
3-»-0.6 


1 7+0.2 
l-'-0.1 


5fi+°- 38 
u - oo -0 .27 


24.23t0-02 


2.6 


1.2/71 




026 




i 9+0.3 
1-0.3 


20.83+3 3* 


1.0/42 


o 7 +0.8 
°-'-0.7 


1 9+0' 2 
x - -0.2 


0-38_ 009 


16.68l0- 9 4 


1.8 


1.0/42 




027 


6 6+ ' 9 


i q+0.2 
1 -°-0.2 


43 22+ 1 ' 18 

-5.46 


0.8/81 


3.8+lj 


1 9+0' 2 


56+ ' 29 
U - OO -0.22 


04 14+1.39 


3.7 


0.8/81 




029 


7q +2.6 
'• y -1.9 


1 s+0-5 

1 - 8 -o.s 


6.55+°;?° 


1.4/25 




1 6+0' 2 


n 90+0.09 
U '^ 8 -0.05 


4 64+ - 30 
^• D4 -1.18 


0.5 


1.3/25 




high 6 


5 4+ - 5 


1 0+ - 1 


25.771^8 


1.0/222 


o 9+O.3 
°' -0.3 


2 1+ 01 
z -!-0.1 


n 07+O.O5 
u -°' -0.04 


21.25+?-^ 


2.3 


1.2/222 




medium 


6 4+°' 6 
-0.5 


1 2+ 01 


16.47«- 2 » 


1.0/212 


o 7+O.4 
°- '-0.3 


2 0+ ' 1 
Z ' U -0.1 


n 09+O.O3 
u ' oz -0.03 


10 i 4 +0.'33 


1.4 


1.1/212 




low c 


6 7 +0 - 9 
D -'-0.8 


i "s+o- 2 


1 09+0. 16 


1.0/115 


3 6+°- 5 


1 7+0.1 
' -0.1 


Q ioq+0.003 
u.lou_ .o 02 


1 oo+0. 04 
l-3»_0.1l 


0.2 


0.8/115 





a Fluxes (corrected for the absorption) are in units of 10 -11 erg cm -2 s — 1 . 
b Luminosities in units of 10 35 erg s — 1 , assuming a distance of 3 kpc. 

c Cash statistics (C-stat) and percentage of 10 4 Monte Carlo realizations that had statistics < C-stat. 
d Blackbody radii are in units of km, assuming the optical counterpart distance of 3 kpc. 

c Intensity-selected spectra. High corresponds to CR > 1 counts s — 1, medium to 0.5 < CR < 1 counts s — and low to CR. < 0.5 counts s~ 

3 RESULTS 
3.1 Light curve 

Fig. [T] shows the 0.2-10keV light curve of IGR J18483— 031 1 of 
the whole campaign after background-subtraction and pile-up, PSF 
losses, and vignetting corrections. Each bin contains a minimum of 
20 source counts. The light curve starts at phase 0.07 , assuming a 
perio d of 18.52 d and an initial epoch MJD 53844.2 dSguera et al.l 
|2007|) . and monitors the flux state through a whole period until the 
following phase 0.51. Superimposed on the long-term orbital mod- 
ulation, variability is seen on short time scales, as shown in the inset 
of Fig.[TJ where a variation by factor of 5.3 in count rate is observed 
to occur in ~ 1.7 hr. This behaviour has been obser ved in several 
SFXTs (e.g.. lSidoli et al.ll2008l ; lRomano et alj2009ah . 

The lowest point in the campaign is a 3a upper limit reached 
on MJD 55006.5-55010.2 at 0.004 counts s" 1 (combined observa- 
tions 017-020, total on-source exposure of 8.3 ks), and corresponds 
to an observed (unabsorbed) flux of 3.9 x 10 -13 (2.1 x 10~ 12 ) 
erg cm ~ 2 s -1 , i f we as sume the XMM-Newton spectrum reported 
by iGiunta et al.l (2009, photon index V = 2.5, absorbing col- 
umn 7Vh = 7.7 x 10 22 cm -2 ). The corresponding luminosity is 
2.3 x 10 33 erg s _1 (assuming the optical counterpart distance of 
3 kpc); to date this is the lowest quiescent X-ray flux value reported 
in the literature for this source. The peak count rate is reached on 
MJD 55016.4 at > 5 counts s _1 ; therefore, the observed dynamical 
range of this source is at least 1200. 

In order to search for spin periodicity, the arrival times of all 
selected events have been converted to the Sol ar System barycentric 
frame, using the BARYCOR code. The Z 2 test iBuccheri et alj 19831) 
on the fundamental harmonics was applied to a sample of source 
photon arrival times for each observation. The search for a timing 
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Figure 2. Spectral parameters as a function of time (see Table fSJ. (a) 
Swift/XKT light curve in the 0.2-10 keV energy band at a day resolution; 
Spectral parameters of the absorbed power-law fit, iVjj, and photon index 
r [(b) and (c)]; Spectral parameters of the absorbed blackbody fit, TVh, 
temperature kT, and blackbody radius [(d), (e), and (f)]. 



feature was performed within the frequency interval 0.01-0. 19 s . 
No presence of coherent pulsations was detected. 
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3.2 Spectra 

Spectra were extracted for each segment in which a detection was 
obtained and a minimum of 120 source counts were available. The 
data were rebinned with at least 20 counts bin -1 to allow x 2 fitting, 
except when the statistics were poor, in which case we adopted 
ICashl l l 19791) statistics and data binned to 1 count bin , instead. 
The simple models we considered were absorbed power laws and 
absorbed blackbodies. The fits were performed in the 0.3-10 keV 
energy band. The results are reported in Table [2] while the spectral 
parameters as a function of time are shown in Fig. [2] In particular, 
for the spectrum of observation 006, the highest in flux during the 
first peak (~ 1.5 x 10 -10 erg cm -2 s _1 ), an absorbed power-law 
model yielded a high absorbing column Nu = (6.3j^'|) x 10 22 
cm~ 2 , and a photon index F — 1.7 ± 0.4 (see Table [2}. As a 
comparison, observation 025, roughly at the same phase, yielded 
N H = (7.lll;o) x 10 22 cm" 2 , and T = 1.6 ± 0.3, hence consis- 
tent with observation 006. Despite the large observed variations in 
flux throughout the campaign, the spectral parameters do not vary 
significantly within the large uncertainties, with the exception of 
the blackbody radii. 

To further investigate the spectral properties of the sources in 
several states, we accumulated all events collected during the cur- 
rent campaign. We extracted events within three intensity levels de- 
pending on count rate, namely, CR > 1 counts s _1 (high, 5447 
counts), 0.5 < CR < 1 counts s _1 (medium, 5264 counts), and 
CR < 0.5 counts s" 1 (low, 2670 counts). We created exposure 
maps for each of these intensity-selected event files and then com- 
bined them (and their exposure maps), and extracted a single spec- 
trum for state. The generation of ancillary response files and spec- 
tral fitting were performed in the same fashion as for the single ob- 
servations. The fit results are reported in Tableland the spectra are 
shown in Fig. [3] Even with the higher statistics afforded by accumu- 
lating all events in three intensity states, no significant variations in 
the column density could be derived. We can confirm, however, that 
the Nn is always in excess of th e Galactic one, 1.4 x 10 22 cm -2 



consistently with that found bv lSguera et alj 1 2007). Similar ly to 
what was found in a sample of 4 SFXTs {Romano et al. 2009^), our 
fits indicate either a hard power law or hot black body. We also note 
that all spectral fits with an absorbed black body resulted in radii 
of the emitting black body region of only a few hundred meters 
(see Tabled, consistent with being emitted from a s mall portion of 
the neu tron star surface, such as its polar caps (see, iRomano et al.1 
l2009ah . 



4 DISCUSSION 

In this paper we report on the first complete monitoring of the X- 
ray activity along an entire orbital period of a Supergiant Fast X- 
ray Transient. This makes these Swift observations a unique data- 
set, which allows us to constrain the different mechanisms pro- 
posed to explain the nature of this new class of X-ray transients. 
IGRJ18483-0311 and IGR Jl 1215-5952 are the only SFXTs 
where both the orbital and spin periods are known, although the 
two systems are very differe nt (P or b ^165 days, P Bp i n ~190 s, in 
IGR Jl 1215-5952; see, e.g.. lRomano et al.ll2009bl) . 

The Swift light curve appears to be highly modulated, with two 
maxima, separated by a time interval consistent with the orbital pe- 
riod of ~18.5 days. A lower limit of 1200 to the dynamical range 
can be obtained from the observed light curve. The different du- 
ration of the two outburst peaks monitored with Swift is probably 




1 2 
Energy (keV) 

Figure 3. Intensity-selected spectroscopy. Upper panel: Swift/XKI data fit 
with an absorbed power law. Lower panel: the residuals of the fit (in units of 
standard deviations). Filled blue circles, green empty circles, and red filled 
triangles mark high, medium, and low state, respectively. 
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Figure 4. Comparison of the Swift/XRJ light curve of IGR J 18483 -031 1 
(crosses) with the prediction of Bondi-Hoyle accretion from a spherically 
symmetric and homogeneous wind. We assumed a distance of 3 kpc. The 
model-dependent orbital phase = corresponds to 54995.83 MJD. 



the result of both a different sampling and a high intrinsic X-ray 
variability. The second peak has a duration of several days, as pre- 
viously observed by INTEGRAL fcguera et alj2007l) . 

The modulation of the overall shape of the light curve with the 
orbital phase can be interpreted as wind accretion along a highly 
eccentric orbit. Thus, we applied different models for the wind ac- 
cretion to gain information on the source parameters. 

The simplest case is a Bondi-Hoyle accretion from a spher- 
ically symmetric and homogeneous wind. Assuming a super- 
iant mass of Mq b = 33 M Q and a radius, Rob ~ 33.8 Rq 
Searle et alJ l2008h . together with the known orbital period of 
Porb = 18.52 d, we tried to account for the overall shape of the 
X-ray light curve, leaving the orbital eccentricity e as a free pa- 
rameter. We assumed the following values for the wind properties: 
a terminal velocity Vao in the range 1400-1800 km s" 1 , a mass 
loss rate M in the rani >e 0.4 -1.5 x 10~ 6 Mq yr" 1 [see Sear le et al.l 



( 2008). lLefever et al.|}2007f) l . In this framework, we obtain the best 
agreement with the observed light curve by assuming an eccentric- 
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prove the agreement between the observed and the calculated light 
curve, we applie d the isotropic clumpy wind model proposed by 
iDucci et alj il2009h. 



The 
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Figure 5. Comparison of the SwiftfXBS light curve of IGR J 18483 -031 1 
(crosses) with the prediction of Bondi-Hoyle accretion from a spheri- 
cally symmetric and homogeneous wind, assuming a magnetic field B = 
4 X 10 11 G. The model-dependent orbital phase <f> = corresponds to 
54997.69 MJD. 



ity e = 0.4, v x = 1800 km s , M = 5 x 10~ 7 M yr _1 , 
(3 — 1. Fig. [4] shows the comparison of the model predictions and 
the observed SwiftfXKT light curve (in units of ergs s _1 , by assum- 
ing a distance of 3 kpc). The model roughly reproduces the shape 
of the X-ray light curve due to the orbital modulation, with the 
largest deviation from the observations being in the time interval 
MJD 55006.5-55010.2, where we observed 4 upper limits. Note, 
however, that we cannot be sure that the low intensity extends for 
4 days continuously, because the four Swift observations consist 
of short snapshots. We investigated the possibility that these up- 
per limits cou ld be due to the onset of a centrifugal inhibition for 
the accretion jDavidson & Ostrikerll 19730 . For the above adopted 
set of wind and orbital parameters, we calculated a new X-ray 
light curve, finding that a low magnetic field of the neutron star 
of B « 4 x 10 11 G could be responsible for the upper limits at 
MJD 55006.5-55010.2 (see Fig.[5j. 

An alternative scenario to explain the low intensity state is an 
X-ray eclipse of the neutron star by the supergiant companion. As- 
suming a circul ar orbit, we obtain a low er limit to the radius of the 
supergiant star teappaport & Josslll983 ) of Rob = 46 Rq , which 
is too large for a B0.5a supergiant JSearle et al.l2008l) . On the other 
hand, our modelling of the X-ray curve suggests a high eccen- 
tricity of at least e = 0.4. Adopting this eccentricity, we derived 
an expected value for the supergiant radius of Rob ~ 39.5 Rq. 
The radii of B0.5Ia stars are usually smaller t han this value, but 
there are several ex ceptions with R &s 40 Rq ( Sear le et alj|2008l ; 
iLefever et alj2007l) . Therefore we cannot exclude that an eclipse is 
responsible for the low luminosity state, in an eccentric orbit. 

Although both centrifugal inhibition and an eclipse can rec- 
oncile the observed low intensity state with the Bondi-Hoyle ac- 
cretion predictions, it is also clear that the spherically symmetric 
and homogeneous wind only reproduces the overall shape of the 
X-ray light curve. It cannot, indeed, account for the very large 
spread around the average behaviour due to the orbital modulation 
and, most of all, the remakable short time scale variability (see, 
for example, the inset in Fig. Q] where a variation by a factor of 
5.3 in count rate is observed in ~ 1.7 hours). The observed short 
time scale variability can be naturally explained by the accretion 
of single clumps composing the donor wind. Thus, in order to im- 
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the effects of accretion from a clumpy wind on the luminosity and 
variability properties of HMXBs. It assumes that a fraction of the 
stellar wind is in the form of clumps with a power law mass distri- 
bution 



(1) 



in the mass range M a -Mb. The rate of clumps produced by the 
supergiant is related to the total mass loss rate Mtot by N c \ = 
^[jy clumps s _1 , where / = M c i/M to t is the fraction of mass 
lost in clumps and < M > is the average clump mass, which can be 
computed from Eq. Q}. Clumps are driven radially outward by ab- 
sorption of UV spectral lines. The following clump velocity profile 
is assumed: w c i(r) = Voa (l — 0.9983^ a )' 3 , where Voa is the 
terminal wind velocity, Rob is the radius of the supergiant, 0.9983 
is a parameter which ensures that v(Rob) ~ 10 km s _1 , and /3 is 
a constant in the range ~0.5-1.5. The model further assumes that 
the clumps are confined by ram pressure of the ambient gas. By 
exploring different distributions for the clump masses and initial 
dimensions, the model can be used to compute the expected X-ray 
light curves in the framework of the Bondi-Hoyle accretion theory, 
modified to take into account the presence of clumps. We sought 
the set of wind parameters yielding the best agreement between the 
calculated and the observed light curve. We found that the observed 
light curve is reproduced well by this wind model by assuming the 
following parameter values: an eccentricity e = 0.4, a mass loss 
rate M to t = 2 x 10" 7 M yr" 1 , = 1800 km s" 1 , /3 = 1, a 
fraction of mass lost in clumps / = 0.75, a mass distribution power 
law index £ = 1.1, a power law index of the initial clump dimen- 
sion distribution 7 = —1 (where Nm c1 oc Rj l clumps s" 1 ), a 
minimum clump mass M a = 10 18 g and a maximum clump mass 
Mb — 5 x 10 2 1 g; moreover, we adop t the force multiplier parame- 
ter obtained bv lShimada et al.1 i 19941) for a B0.5Ia star (k = 0.375, 
a = 0.522, S = 0.099). 

Fig. [6] shows the comparison of the Swift/XKT light curve of 
IGR J18483— 031 1 with the isotropic clumpy wind model predic- 
tion. Further acceptable solutions can be found by assuming wind 
parameters £, /, 7 in the allowed ranges plotted in Fig. [7] and e = 
0.4 ±0.1, Mtot = (2±1) x 10- 7 M Q yr-Suoo = 1800 km s" 1 , 
/3 = 1, 10 18 < M c i ^ 10 21 g. As Fig.[7]demonstrates, the com- 
parison of the observed light curve with the clumpy wind model 
allowed us to constrain the parameters responsible for the degree 
of inhomogeneity of the wind. In particular, we found that a very 
large fraction of the mass lost from the supergiant is contained in 
the clumps (0.7 < / < 0.78), and we obtain the value of C, (which 
controls the shape of the clump formation rate distribution) with an 
accuracy of ^15%. 

Fig. [6] shows that the peak luminosities, the dynamic range 
involved by the flares, and the orbital modulation and the low lu- 
minosity state (MJD 55006.5-55010.2) observed are reproduced 
well by the clumpy wind model, even without invoking either a 
centrifugal barrier or an X-ray eclipse. Indeed, from the calcu- 
lated light curve, we determined that the probability to observe the 
source at the inter-clump luminosity level in the range of phase 
0.2 < <j> < 0.8 is ~ 25%: with the binomial distribution func- 
tion, we obtain a probability to measure 4 low luminosity states 
of ~ 5%. Therefore, the upper limits can be explained with the 
accretion of the intra-clump wind with a low density, even with- 
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Figure 6. Comparison of the Swift/XKT light curve of IGR J 18483 -03 11 
(top) with the pr ediction (bottom) of the new clumpy wind model of 
iDucci et al] (2009). The model-dependent orbital phase <f> = corresponds 
to 54995.83 MJD. 



out invoking centrifugal inhibition or an eclipse. Finally, we note 
that the win d parameters we o btain applying our spherical clumpy 
wind model jDucci et alj2009h are very similar to those explaining 
the Vela X-l X-ray light curve. Indeed, the two systems have very 
similar donor stars. 
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Figure 7. Graphs of the allowed parameters f , /, 7 (filled region), obtained from the comparison between the observed and the calculated light curves. 



